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The existed interpolation method, based on the second-order tetration 
polynomial, has the asymmetric property. The interpolation results, for each 
considering region, give individual characteristics. Although the 
interpolation performance has been better than the conventional methods, 
the symmetric property for signal interpolation is also necessary. In this 


paper, we propose the symmetric interpolation formulas derived from the 
second-order tetration polynomial. The combination of the forward and 
Keywords: backward operations was employed to construct two types of the symmetric 
interpolation. Several resolutions of the fundamental signals were used to 
evaluate the signal reconstruction performance. The results show that the 
proposed interpolations can be used to reconstruct the fundamental signal 
Hyperpower and its peak signal to noise ratio (PSNR) is superior to the conventional 
Signal reconstruction interpolation methods, except the cubic spline interpolation for the sine wave 
Super exponentiation signal. However, the visual results show that it has a small difference. 
Moreover, our proposed interpolations converge to the steady-state faster 
than the cubic spline interpolation. In addition, the option number increasing 
will reinforce their sensitivity. 
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1. INTRODUCTION 

Recently, the technology based on digital information is rapidly growing up. It can change the 
human lifestyle. Some information can be easily learned and applied to the real world. However, all data 
cannot be kept in digital form due to a lot of limitations, such as storage size, sampling technique, or 
communication constraints. In order to obtain the unknown data from its related known data, interpolation is 
an important method. It has been widely applied in several applications, i.e., signal reconstruction [1]-[9], 
signal and image denoising [10]-[12], medical image processing [13]-[18], image resizing [19]—[28], image 
classification [29], image compression [30], image enhancement [31], change detection [32], video 
processing [33], and more [34], [35]. The interpolation is used to generate the new data points from the 
discrete set of known data points in case of its position is within the range of known data. The simplest 
method is the nearest neighbor interpolation employed the data value of the nearest location to estimate the 
new data points. This method is the same fashions of the piecewise constant function. It is an advantage for 
the lowest computing time but lacks performance for other situations. These problems can be solved by 
substitute the piecewise constant function with the higher-order polynomial, i.e., the linear interpolation [36] 
or the cubic interpolation. The linear interpolation gives the computing time lower than the cubic 
interpolation, but the reconstruction performance in smoothing signal is poor. We tried to apply the tetration 
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mathematics, which is the next mathematical operation after exponentiation in our previous work, for 
modeling the alternative formulas for signal interpolation [2] as: 


Tetration polynomial: the tetration is one of the mathematical operations located after exponentiation but 
before pentation. It was firstly proposed by Goodstein [37] in 1947. He used the notation "x instead of 
the application n — 1 time of x, that is: 


n 


oo 
x 


ny = a” > (1) 


however, there are many notations for the tetration, i.e., x TT n for [38] and uxp% for [39]. Moreover, it 
also can be called for a different name such as superexponentiation [40], hyperpower [41]. This 
mathematical operator was applied in the problem of compacting Church numerals [42], and our previous 
work [2]. The form of the traditional polynomial function was used to model the tetration polynomial 
function. It is the additional series that each term is the multiplication of the arbitrary constant (ap) and 
any order tetration. We can write in the summation form as (2): 


f(x) = n=0 an (°x), m E N (2) 


Existed tetration interpolation: the simple quadratic tetration polynomial was considered in our previous 
work [2], written as (3): 


fr(X) =a +a x +a: *x (3) 


If we know the data of any considering signal as (x1, Y1), (X2, Y2), and (x3,y3), therefore we can 
construct the system of equation for these points as (4): 


2 
1 x, 


xX, Ap 
1 X2 2%, e 
1 x3 xq) 1% 
In the tasks of discrete signal processing, we can define x,,xz and x3 as the non-negative integer 
numbers. Thus, the solution of above system of equation can be found. Several techniques can be used for 


calculating the parameters dj, a4, and az. In this paper, we employed the Cramer's rule for determining 
solution of (4) as (5). 


yı 
y2 


Y3 


(4) 
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( 2xgX2- 2x9x3)+( 2x x3- 2x3X1)+( 2x1 - 2x4Xx2) 


Forward and backward operations: The above interpolation formula (3) has the shift variant and 
asymmetric properties [2]. It gives an individual value for each domain range. For example, if we set the 
value of y,, yz, and y3 as any fixed real number and set x1, X2, X3 are the non-negative integer number 
with arranged order. The interpolation results can be shown as Figure 1(a). The word Option means as the 
set of sampling step position. Each option can be defined as x, = 0,x2 = 1, x3 = 2 for Option 1; x, = 1, 
X, = 2, X3 = 3 for option 2; x, = 2, x} = 3,x3 = 4 for option 3, and so on. The interpolation results 
were compared with the curve fitting of the quadratic polynomial. The simulation results show that the 
interpolated curve of the quadratic tetration polynomial for each option changes with asymmetry. It gives 
a high sensitivity on the right and provides a low sensitivity on the left. The effect of the option is for 
reinforcement of the response varied on the number of the option. The formula (3) will be called forward 
operation because it was operated from left to right of considered points. In the other hand, if the 
positions are operated from right to left, we called backward operation, which the interpolation results can 
be shown as Figure 1(b), and can be expressed as (6): 


fp(X) = ao + a ' (xo — x) + az: 2 (xo — x) (6) 
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where ap, @,, and a, are the parameters, which can be calculated from (5), and xo is a backward constant 
being equal two times of the option number. 

Following the above methods, we found that it can be applied to reconstruct the fundamental 
degraded signals with many specific properties, i.e., shift variant and asymmetric properties. However, the 
symmetric property is also necessary for several applications. In this paper, the symmetric signal 
interpolation based on a quadratic tetration polynomial was proposed. The rest of this paper can be organized 
as follows. The methods for constructing the interpolation formulas, both asymmetric and symmetric 
formulas, are presented and described in section 2. Then, the research method is illustrated in section 3. The 
experimental results of the fundamental signal are demonstrated and discussed in section 4. Finally, we give 
the concluding remarks in section 5. 
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Figure 1. Curve fitting of the quadratic tetration polynomial of option 1, 2, and 3, compared with 
The normal quadratic polynomial for (a) forward operation and (b) backward operation 


2. PROPOSED METHOD 
2.1. Asymmetric tetration interpolation 

The option 1 was selected to demonstrate our approach techniques. That is, the value of x1, x2, and 
x3 are 0, 1 and 2, respectively. In this option, the (5) can be reduced in the short form as 


ao = = (2% + 2y¥2—y3), a, = “(y +yz), and a, = o — 2y, + y3). Therefore, we can build the 
Forward interpolation formula by substitute of the constant ao, a, and az into (3). That is: 


Yrorward = A1Y1 + Azy2 + A3y¥3 (7) 


where 
1 2 
A; =3(2-3x+ x) 


A2 


+(2+3x-2(7x)) ,x € (0,2) 


A3 


1/_ 2 
Likewise, the backward interpolation formula can be expressed as (8): 


Yeackward = B1Y1 + Boyz + B3Y3 (8) 


where 


B, = z (-1+ ?(xo — x)) 
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B, =1(2 +3- x)-2: 2(xo —x)) , x € (0,2) and xọ = 2. 


B; = E2 — 3(xo — x) + ?(xo — x)) 


The interpolation results on the left side (0, 1] and the right side [1, 2) give the different results 
(asymmetric property). We can demonstrate four different techniques of asymmetric interpolation as follows. 
Firstly, the interpolation result on the left side of the Forward operation is considered. It was called forward- 
left (F-L) tetration interpolation. Also, if the interpolation result on the right side of the forward operation is 
considered, it will be called forward-right (F-R) tetration interpolation. Some behavior of such techniques for 
square wave signal reconstruction can be demonstrated in Figure 2(a) and Figure 2(b). In the same way, the 
backward-left (B-L) tetration interpolation and the backward-right (B-R) tetration interpolation will employ 
the interpolation result of the backward operation on the left and right sides, respectively. Figure 3(a) and 
Figure 3(b) illustrates the reconstruction of the B-L and B-R tetration interpolations for the square wave 
signal case. 

We can summarize the preliminary dominant results as follows. The F-R and B-L interpolation 
techniques give a higher edge overshoot than the F-L and B-R techniques. The edge overshoot does not 
appear in both terminal sides of the changing region that is different for the cubic spline interpolation. The 
edge overshoot of the same operation does not appear on the same location, but it gives a similar location 
when they have the same considered side. For example, the F-L and B-L are the same edge overshoot 
location. Likewise, the F-R and B-R have the same location of the edge overshoot. 


2.2. Symmetric tetration interpolation 

In order to correct the asymmetric behavior, the arithmetic mean of the combination of the F-L 
interpolation and B-R interpolation is used to formulate the symmetric interpolation formulas. Because of the 
behavior of both F-L interpolation and B-R interpolation, it is low sensitivity with the high-frequency 
components or the changing region as shown in Figure 2(a) and Figure 3(b). We will call this interpolation as 
low sensitivity tetration (LOST) interpolation. This formula can be simply expressed as follows, 
LOST interpolation; 


Yost = 4 Vi-1 F A2Yi + A3Yi+1 + A4Yi+2 (9) 
where 


a, = Gi + ?(1-x)) 


1 

ar=z(7-6:x+ 2x — 2: ?(1— x)),x € (0,1) 
1 2 2 

az =z(1+6:x-2: x+ 7(1-x)) 


1 
a, ==(-1+ °x) 
6 
Likewise, the arithmetic mean of the combination of the F-R interpolation and B-L interpolation is 
used to formulate the high-sensitivity symmetric interpolation formula. This formula has the opposite 
behavior with the LOST interpolation because it was constructed from two high sensitivity formulas as 


shown in Figure 2(b) and Figure 3(a). Therefore, this formula will be called high sensitivity tetration (HIST) 
interpolation. Its formula can be expressed as follows. HIST interpolation; 


Yuist = Piyi-1 + Boyi + P3Yi+1 + BaYire (10) 


where 


B= 5(-A-3-x4 (1 + x)) 
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Be = (4+3 x-2: ?(1 +x) + ?(2 -—x)),x € (0,1) 
1 
b= VOA Arn 72-2) 


Ba =E(-4+3:x4 2(2 — x)) 
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Figure 2. Comparison of the square wave interpolation result with (a) F-L tetration interpolation and 
(b) F-R tetration interpolation 
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Figure 3. Comparison of the square wave interpolation result with (a) B-L tetration interpolation and 
(b) B-R tetration interpolation 


3. RESEARCH METHOD 

In order to evaluate our proposed method, the MATLAB program version R2018a was employed to 
implement the experimental results. Three fundamental synthetic signals i.e., square wave signal, a saw-tooth 
wave signal, and sine wave signal with several sampling resolutions were simulated and applied. It was 
designed for observing the performance limitations. Two comparison types were demonstrated as the quality 
metric and visual results. For the metric results, the peak signal to noise ratio or (PSNR) was employed to 
indicate the reconstruction quality. Our proposed methods were compared with three conventional methods, 
i.e., nearest-neighbor interpolation, linear interpolation, and cubic-spline interpolation. For the vision results, 
both our methods were illustrated in graph plot form and were compared with the cubic spline interpolation. 
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4. SIMULATION RESULTS AND DISCUSSION 

The experimental results were devided into two parts as follows. The fundamental synthetic signals 
were employed to test our method performance in the first part. After that, the behavior of the option was 
demonstrated for observing the trend of the experimental result in second part. 


4.1. Fundamental synthetic signal testing 

The synthetic square wave signals were generated with 5, 15, 25, 35, and 45 samples per period and 
then were applied to all comparison methods, including our methods. The quality metric results were 
illustrated in Table 1. We can see, the LOST interpolation gives the highest PSNR results for all resolutions. 
For the HIST interpolation, it gives the PSNR value lower than others, except the nearest-neighbor method. 
The vision results were shown in Figure 4. Our methods were plotted comparing with the original data and 
the cubic spline interpolation. The HIST interpolation gives the behavior similar to the cubic spline 
interpolation for the edge region, but it tilts out of the edge as shown in Figure 4(a) in the zoom image. 

Moreover, the HIST gives a behavior similar to the linear interpolation when it is out of the edge 
area. This behavior may be from the tetration polynomial that also has the first order normal polynomial 
(linear term) and its combination techniques. The vision results of LOST interpolation, it was demonstrated 
in Figure 4(b). We can see, the edge ringing is lower than others. Moreover, it also has a similar property 
with the HIST interpolation of both the tilt characteristic and off-edge region behavior. 


Table 1. The PSNR value versus the sampling resolutions for the square wave signal 
Resolution (Samples/Period) _ Nearest _ Linear Cubic LOST HIST 
5 13.9 15.2 14.7 15.2 15.0 


15 18.6 20.1 20.1 20.2 20.0 
25 20.7 22.3 22.3 22.4 22.2 
35 22.3 23.8 23.8 23.9 23.6 
45 23.4 24.9 24.9 25.0 24.7 
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Figure 4. Comparison results between the cubic spline interpolation, (a) HIST interpolationand and (b) LOST 
interpolation for the square wave signal 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1893-1903 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1899 


In the second experiment, we apply the saw-tooth wave signal to test our method performance. The 
visualization results of its behavior are demonstrated in Figures 5(a) and 5(b). For the quantitative 
comparison, the PSNR value of each method were shown in Table 2. The PSNR values of the LOST 
interpolation are also higher than others, following with the HIST interpolation. The zoom image in 
Figure 5(a) indicates that the HIST interpolation is also similar to the cubic spline interpolation, but it tilts out 
of the instantaneously changed edge. For this behavior, we may summarize as the tilt direction is opposite 
with the edge that has a higher slope. Moreover, the results also show that the linear property gives the 
converges to the steady-state faster than the cubic spline. 
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Figure 5. Comparison results between the cubic spline interpolation, (a) HIST interpolation 
and (b) LOST interpolation for the saw-tooth wave signal 


Table 2. The PSNR values versus the sampling resolutions for the saw-tooth wave signal 
Resolution (Samples/Period) Nearest Linear Cubic LOST HIST 


5 14.3 16.3 16.4 16.6 16.5, 
15 18.9 21.2 21.3 21.5 21.4 
25 21.1 23.4 23.5 23.7 23.6 
35 22.5 24.9 25.0 25.2 25.1 
45 23.6 26.0 26.1 26.3 26.2 


For the third experiment, we applied our interpolation methods to reconstruct the sine wave signal. 
The interpolated results can be illustrated in Table 3 for the quality metric results and in Figures 6(a) and 6(b) 
for the visual results. The results show that the cubic spline interpolation gives a higher PSNR than others, 
following by the LOST and HIST interpolation, respectively. For vision result of the HIST interpolation, it 
confirms that the tilt direction is dependent on the slope between both sides of the interpolation interval. The 
HIST interpolation is more sensitive with a rate of data change than the LOST interpolation. Although the 
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quality metric of the LOST interpolation looks more different than the cubic spline interpolation, the visual 
results are too close to cubic spline interpolation. 


Table 3. The PSNR values versus the sampling resolutions for the sine wave signal 


Resolution (Samples/Period) Nearest Linear Cubic LOST HIST 
5 21.9 30.2 48.2 37.9 38.8 
15 31.4 49.0 92.6 62.3 51.7 
25 35.8 57.8 114.0 72.1 60.1 
35 38.7 63.7 127.1 78.2 65.9 
45 40.9 68.1 136.4 82.7 70.2 
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Figure 6. Comparison results between the cubic spline interpolation, (a) HIST interpolation and (b) LOST 


4.2. Option behavior for symmetric tetration interpolation 


increasing will provides higher sensitivity than the lower value. 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1893-1903 


interpolation for the sine wave signal 


In this subsection, the option 1, 2, and 3 were demonstrated compared with the original square wave 
signal and the result of the cubic spline interpolation for observing the behavior of our symmetric tetration 
interpolation. The vision results, both LOST and HIST interpolations, were illustrated in Figures 7(a) and 
7(b). The results show that the edge overshoot of the LOST interpolation decrease when the option number is 
higher. Vice versa, the edge overshoot of HIST interpolation increases following the option number. For 
observing the quality metric with several sampling resolutions, the graph plot of the PSNR of the LOST, 
HIST, and cubic spline was presented in Figure 8. The LOST interpolation also gives a higher PSNR than 
others. The higher option number gives a better performance than the lower option number. However, the 
trend of performance result is a more similar value when the option number is much higher. Anywise, the 
higher option number gives the poor performance for the HIST interpolation. That means the option number 
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Figure 7. Comparison results between the cubic spline interpolation, (a) HIST interpolation and (b) LOST 
interpolation for the square wave signal with options 1, 2, and 3 
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Figure 8. The PSNR comparison of HIST and LOST interpolations with options 1, 2, and 3 compared with 
the cubic spline interpolation for the square wave signal 


5. CONCLUSION 

In this paper, the symmetric tetration interpolation methods were proposed. It derived from the 
second-order (quadratic) tetration polynomial. Four types of asymmetric interpolation formulas were 
summarized and investigated. Then, they were analyzed, demonstrated, and discussed for testing with square 
wave signal only. After that, the symmetric formulas were approached into two types, such as the HIST 
interpolation and the LOST interpolation. The combination of the Forward and Backward operations was 
used to formulate the above methods. The typical synthetic signals, i.e., square wave, saw-tooth wave, and 
sine wave signals, were applied to evaluate the interpolation performance. Three conventional interpolation 
methods, i.e., nearest neighbor interpolation, linear interpolation, and cubic spline interpolation were 
employed for comparing the proposed method performance with the PSNR indicator. The simulation results 
show that the HIST interpolation has a high sensitivity to the data change. In contrast, the LOST interpolation 
gives a low sensitivity to the data change. The LOST interpolation gives the PSNR value, for all tested 
resolutions, higher than others, except for the sine wave signal that is lower than the cubic spline 
interpolation. However, the vision results show that it has a small difference. The option number increasing 
provides higher sensitivity than the lower value. Moreover, both our proposed symmetric methods also give 
the tilt behavior, which is dependent on the slope of sample data for both sides of the interpolation interval. 
Finally, both our symmetric interpolations still converge to the steady-state faster than the cubic spline 
interpolation. 
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